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Abstract. The cumulative light curves of a large sample of gamma ray bursts (GRBs) were obtained by summing 
the BATSE counts. The smoothed profiles are much simpler than the complex and erratic running light curves 
that are normally used. For most GRBs the slope of the cumulative light curve (S) is approximately constant 
over a large fraction of the burst. The bursts are modelled as relaxation systems that continuously accumulate 
energy in the reservoir and discontinuously release it. The slope is a measure of the cumulative power output of 
the central engine. A plot of S versus peak fiux in 64 ms (P64ms) shows a very good correlation over a wide range 
for both long and short GRBs. No relationship was found between S and the GRBs with known redshift. The 
standard slope (S'), which is representative of the power output per unit time, is correlated separately with P64ms 
for both sub-classes indicating more powerful outbursts for the short GRBs. S' is also anticorrelated with GRB 
duration. These results imply that GRBs are powered by accretion into a black hole. 
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1. Introduction 

Cosmological gamma ray bursts (GRBs) em it an extraor 



dinary amount of energy in gamma rays ( Costa et al 
|1997| ; Ivan Paradijs et al. , 1997| ; |Piran, 1999|) . The source of 



this energy may be a cataclysmic event involving mergers 
of compact objects such as neutr on star binaries or neu- 
tron stars and black hole binaries ( RufFert fc Janka, 1999 ) 
or the formation of a black hole during or after the col- 



2. Timing analysis of the GRB profiles 

A large sample of 498 of the brightest BATSE GRBs with 
data combined from the four energy channels was anal- 
ysed and denoised using either wavelets or median filters. 
A detailed account of this process including sample selec- 
t ion and pulse analysis in GRBs has been give n elsewhere 
( iQuiUigan et al., 2002| ; [McBreen et al., 200l|) . The GRB 



lapse of massive stars (Rees fc Meszaros, 1994 ; Paczynski, 



1998; Vietri fc Stella, 1998; MacFadyen fc Woosley, 1999; 



Reeves et al., 2002 ). GRB light cu rves are complex and er 



ratic ( Fishman fc Mcegan, 1995 ). There are a number of 
recent results on the properties of the pulses in GRBs and 
their relationship to the duration Tgo (Ramirez-Ruiz & 
Fenimore, 2000; Salmonson, 2000; Norris, 2002; McBreen 
et al., 2002a; Gupta et al., 2002; Quilligan et al., 2002 and 
references therein). The correlated pulses in GRBs have 
a unique set of properties. In this paper we show in sec- 
tions 3 and 4 that the slope of the cumulative light curve 
of most short and long GRBs is approximately linear and 
the bursts can be modelled as relaxation systems. In ad- 
dition the slopes are highly correlated with the peak flux 
and anticorrelated with GRB duration. The large and uni- 
form BATSE sample of GRBs were used in this analysis 



(Fishman fc Meegan, 1995; Kouveliotou et al., 1993) 
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sample was restricted here to bursts that had the on-board 
summed count from only two Large Area Detectors and 
includes 55 GRBs with Tgo < 2 s analysed at 5 ms res- 
olution, 250 GRBs with Tgo > 2 s at 64 ms resolution. 
To extend the sample to include additional GRBs with 
Tgo > 100 s, a further 71 GRBs were included and anal- 
ysed at 256 ms resolution. The cumulative light curve was 
obtained by converting the BATSE rates to counts and 
taking the cumulative sum. The cumulative count from 
the same number of on-board detectors is uncorrected 
for the detector response and the angles between the de- 
tectors and GRB sources. The errors from these effects 
should be less than ±30%. As a check on this procedure, 
the cumulative count for the GRBs was compared with 
the BATSE fluences. There is a high degree of correlation 
between the two quantities showing that the cumulative 
count from the same number of detectors is a good mea- 
sure of the burst. There is close agreement with the cumu- 
lative light curves o btained in the spectroscopic sample of 
Preece et al. (2000|) . 
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3. Results 

The running and cumulative light curves of a sample of 
bursts, that cover a wide range in T90, are given in Fig. 1. 
The cumulative light curves smooth the running profiles. 
The slope of the cumulative light curve (S) was measured 
for the section of the GRB that included more than 50% 
of the total counts and could be represented by a straight 
line, with a precision in the measured slope that was al- 
ways much better than a factor of 2. Not all GRBs were 
fit with a single straight line because of long time inter- 
vals with very little or no emission (Fig. Id) and some 



of these GRBs have been discussed by Ramirez-Ruiz & 



Merloni (2001). In these cases two or more straight lines 
were used to fit the separate periods of emission. The av- 
erage value of the slopes was adopted for that GRB. This 
procedure was used for 66 GRBs with Tgg > 2 s and 2 
GRBs with T90 < 2 s. The cumulative light curves of a 
number of GRBs have measurable curvature and in these 
cases a curve is a more accurate fit than the straight line 
adopted here. These GRBs are the subject of a separate 



pubhcation (McBreen et al., 2002b) 



The median percentage of the integrated counts over 
which the slope was measured is 71% for GRBs with T90 > 
2 s and 83% for short GRBs. The corresponding values of 
the median percentages of T90/T50 are 34%/127% and 
68%/140%. T50 is a better measure than T90 of the active 
period of the burst over which the slopes were measured. It 
is also better anticorrelated with the slope than T90 (Table 
1, Fig. 2c). This effect is illustrated in Fig. Id where 75% 
of the integrated counts is used for the slope over only 
24% of T90 whereas in Fig. la the corresponding values 
are 69% and 59% respectively. There is a long quiescent 
interval in the GRB in Fig. Id which contributes to the 
small percentage value of T90. 

The slopes of the cumulative light curves of 376 GRBs 
are plotted in Fig. 2a versus the peak fiux in 64 ms 
(P64ms)- There is a good correlation between S and P64ms- 
Spearman rank order correlation coefficients p, associated 
probabilities and index of best fit power law for these 
quantities are listed in Table 1. The standard slope S' was 
obtained by dividing the slope of the cumulative counts 
(S) by the time over which the slope was measured and 
hence it represents the cumulative counts or power output 
of the source per unit time. The standard slope is plotted 
versus P64ms in Fig. 2b. The correlation coefficients are 
listed in Table 1 for short and long GRBs. The best fit 
power laws have indices close to 1.6 for both classes of 
GRBs but there is considerable spread in the values with 
the short bursts displaced from the long bursts by ~ 10^. 

The standard slope is plotted versus T50 in Fig. 2c. 
The correlation coefficient is -0.90 with a best fit power 
law index of -1.1. The corresponding values for T90 are 
-0.85 and -1.08 respectively. The values are also listed in 
Table 1 for T90 and T50 versus S. S' is better anticorrelated 
with duration than S. 

Fenimore (1999) obtained the average temporal pro- 
file of 98 GRBs by normalising to a standard duration 



Table 1. Spearman rank order correlation coefficients, p, 
the associated probabilities and index of the best fit power 
law for a range of GRB properties. In all cases the high 
values of p show either strong correlations or anticorrela- 
tions between the burst properties. 



Properties 



Probability 



Power Law 
Index 



P64ms VS. S 


0.84 


< io-^« 


1.04 


P64ms VS. S' (T90 > 2 S) 


0.72 


< io-*« 


1.63 


P64ms VS. S' (T90 < 2 S) 


0.69 


6.1 X 10"^ 


1.65 


T90 VS. S' 


-0.85 


< io-^« 


-1.08 


T90 VS. S (T90 > 2 s) 


-0.58 


1.7 X 10"^" 


-0.54 


T90 vs. S (T90 < 2 s) 


-0.39 


3.0 X lO-^* 


-0.37 


T50 vs. S' 


-0.9 


< io-*« 


-1.10 


T50 vs. S (T90 > 2 s) 


-0.64 


1.9 X 10"^^ 


-0.53 


T50 vs. S (T90 < 2 s) 


-0.53 


3.7 X 10"^ 


-0.32 


Fluence vs. 


0.92 


< io-^« 


0.92 


Cumulative counts 



and a standard peak counts. The average profile is rea- 
sonably fiat over more than 50% of the duration (Fig. lb 
in Fenimore (1999)) implying a linear increase in the cu- 
mulative profile and hence in agreement with the results 
presented here. 

4. Discussion 

4.1. GRBs as relaxation systems 

It can be assumed that the sum of the counts in the bursts 
(Fig. 1.) is a good measure of the integrated energy emit- 
ted by the source because the peak energy li es well within 
the BATSE band ( |Fishman fc Meegan, 199^ ). Many mod- 
els of GRBs consist of a newly formed black hole that ac- 
cretes fro m a remnant torus that is cooled by neu t rino 



emission (Popham et al., 1999 ; Narayan et al., 2001 ; Lee 



|fc Ramirez-Ruiz, 2002| ). The energy to drive the relativis- 
tic jets and bursts may be extracted from the disk and 
spinning black hole by MHD processes and neutrino anni- 
hilation. The efficiency of these processes in creating rela- 
tivistic jets is on the order of one percent (MacFadyen & 
Woosley, 1999| ) 

Models of this type can be usefully compared with a 
relaxation system ( Palmer, 1999| ) which is taken to be one 
that continuously accumulates energy from the accretion 
process and discontinuously releases it. The energy in the 
reservoir at any time t is 



E{t) ^Eo+ I R{t)dt - ES*, 

J O 



(1) 



where Eq is the energy stored in the reservoir that accu- 
mulates energy at a rate R(t) and discontinuously releases 
events of size Si. 

The simplest system is referred to as a relaxation os- 
cillator where there is a fixed level or trip-point that trig- 
gers a release of the energy when E — Emax. The sof t 
^amma ray repeater (SGR) (|Palmer, 199S|; [Gogiis et al 
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2000; Hurley et al., 1994) and GRB pulses are not con- 
sistent with this oscillator. More complicated behaviour 
occurs when the accumulation rate, trigger rate or release 
strength are not constant. If the system starts from a min- 
imum level E = Emin, accumulates energy at a constant 
rate R = r, the sum of the releases is approximately a lin- 
ear function of time i.e. HSi oc rt. This model can account 
for the approximately linear increase in cumulative counts 
from GRBs (Fig 1). The pulses in GRBs have a tendency 
to keep the cumulative count close to a linear function and 
maint ai in a steady state situation , 



i: 



3- 



Ramirez-Ruiz & Merloni (2001) found a correlation 
between the duration of an emission episode in a multi- 
peaked burst and the duration of the preceding quiescent 
time which is similar to the above scenario. The system 
could build up its energy probably via an MHD instabil- 
ity driven dynamo and reach a near critical or metastable 
level. A local instability could cause a rapid dissipation 
of all the stored energy. The system will tend to return 
to a more stable configuration characterised by a certain 
threshold energy Eq, or a sub-critical magnetic field config- 
uration. The source then becomes quiescent. Interestingly 
both models can be unified by noting that each time the 
system is completely drained by a total release of accumu- 
lated energy, the central engine goes quiescent but other- 
wise the energy extraction is usually in episodes that are 
incomplete releases of energy. In this case the longer the 
quiescent time, the higher the stored energy from the next 
episode. Such a situation may give rise to the observed cor- 



relations between long quiescent times (Ramircz-Ruiz & 
Merloni, 2001), correlated pulse properties and intervals 
betwe e n pulses (iQuilligan et al., 2002 ; McBreen et al. 
2002a| ; [Nakar fc Piran, 2002]) . This is a different mecha- 



nism from any relaxation oscillator which forces a release 
of energy when the system reaches an upper level. 



4.2. Relationships between the slopes and peak flux 

There is a significant correlation, that extends over a range 
of ~ 10^, between P64ms and S (Fig. 2a and Table 1). The 
seven GRBs, detected by BATSE, with measured redshifts 
are also plotted in Fig. 2. These sources cover a wide range 
with no obvious relationship between S, P64ms and z. 

The separation of GRBs into two sub-classes was char- 



acterised by durations > 2 s and < 2 s (Kouveliotou et al.. 
1993| ). The plot of the standard slope versus P64ms (Fig. 



2b) also separates the GRBs into two classes with the short 
GRBs having a more powerful output power per unit time. 
The cosmological GRBs with known redshift have values 
of S and S' that range from 7 x 10'*'' to 5 x lO^^ Watts 
and 1.5 x 10'''' to 2 x 10''^ Watts per second respectively. 
There is no evidence for a new class of GRBs with a dif- 
ferent relationship between the standard slope and P64ms- 



4.3. Relationships between the slopes and durations 

The data presented in Fig. 2c shows that as the stan- 
dard slope increases the values of T50 decrease and this 
effect is present in both sub-classes (Table 1). The trend is 
quite revealing and shows that the smaller the value of T50 
the greater the standard slope or the cumulative power 
output per second from the source. The standard slope 
plays an important role in determining T50. The median 
values of the pulse properties and time intervals between 
pulses were found to increase with Tgo (McBrecn ct al 



2002a). The opposite relationship exists here between S' 
and burst duration (Table 1) implying that as the stan- 
dard slope increases there is a corresponding decrease in 
the pulse properties. GRBs with high accretion rates have 
large values of the standard slope, fast pulses and short 
durations whereas lower accretion gives lower values of S', 
slower pulses that are further apart and larger values of 
T90/T50. These results provide strong evidence that GRBs 
are powered by hyperaccretion into a black hole from a 



standard type engine ( Salmonson, 2GGC ; 


Frail et al., 2001 


Panaitescu & Kumar, 2001 




Piran ct al., 2001 


)■ 
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Fig. 1. The running (dashed) and cumulative (solid) light 
curves of the BATSE bursts with trigger numbers a) 3128, 

b) 3057, c) 3042, d) 7560 , c) 2217 with count per 64 ms 
and cumulative count scales on the left and right vertical 
axes. The insert gives the straight line fit (dashed) to the 
cumulative count (solid) for the relevant scction(s) of the 
GRB. The vertical axes in the inserts are the normalised 
cumulative count. 
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Fig. 2. The values of P64ms are plotted versus a) the slope 
S and b) the standard slope S' of the GRB cumulative hght 
curve for three categories of GRBs i.e. Tgo < 2 s (red), 
Tgo > 2 s (green) and the additional sample with Tgo > 
100 s (blue). T50 is plotted versus S' in c) for the same 
three catagories. The seven GRBs with known redshift and 
detected by BATSE are labelled (van Paradijs et al., 2000; 
Castro- Tirado, 2001 and references therein). The BATSE 
trigger numbers and redshifts are given in the top figure. 
An extension of the peak fiux limited sample with T90 > 
2 s to lower values should populate the region containing 
GRBs 1 and 2 with known z in b). 



